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Abstract 

SuperB experiment could represent an ideal environment to test a 
new U{1) symmetry related to light dark forces candidates. A promis- 
ing discovery channel is represented by the resonant production of a bo- 
son [/, followed by its decay into lepton pairs. Beyond approximations 
adopted in the literature, an exact tree level calculation of the radia- 
tive processes e^e~ ^ 7, [/ ^ /i^/i~7, e^e~7 and corresponding QED 
backgrounds is performed, including also the most important higher-order 
corrections. The calculation is implemented in a release of the generator 
BabaYaga@NLO useful for data analysis and interpretation. The distinct 
features of U boson production are shown and the statistical significance 
is analysed. 

1 Introduction and theoretical framework 

The aim of this paper, which is based on [1], is to illustrate Dark Matter models 
which imply the existence of a new vector boson, which carries a new dark force 
and which is lightly coupled to the photon. In particular a Monte Carlo event 
generator useful to describe this kind of light dark forces signatures at leptonic 
colliders will be presented and potentiality of a SuperB experiment for this kind 
of studies will be shown. 

The existence of an abelian gauge simmetry with an associated light U boson 
which can interact with a really small coupling with Standard Model (SM) 
particles has been proposed by a wide class of new physics models[ - ]. 

An important phenomenological support to this class of models comes in re- 
cent years, when standard astrophysics and particle physics seem fail to explain 
striking experimental signals in terms of known sources. 

These signals are represented by the 511 KeV gamma-ray signal from the 
galactic center observed by the INTEGRAL satellite[ ], the excess in the cosmic 
ray positrons reported by PAMELA [ ]. the total electron and positron flux 
measured by ATIC[ ], Fermi[ ], and IIESS[ ], the annual modulation of the 
DAM A/LIBRA signal [ ] and the features of the low-energy spectrum of rare 
events reported by the CoGeNT collaboration[13, 14]. 

If a new secluded gauge sector under which the SM particles remain un- 
charged is included into the theoretical description, these evidences can be 
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Figure 1: Typical reaction useful at leptonic colliders to study U boson proper- 
ties. 

comprehensively interpreted in terms of WIMP Dark Matter (DM) particles 
connected to SM ones with interaction terms varying from model to model. 

The simplest assumption [ ] is to add an extra U{1) symmetry to SM sym- 
metry group which describe a new dark force and suppose this force as carried 
by a new vector boson (U). It is possible to suppose also that U boson can 
communicate with the SM through a kinetic mixing term of the form 



describing the interaction of the U boson with SM photon. This mixing term 
could occour through loop effects due to really massive WIMPs both coupled 
to ordinary photon and U boson. In this case the e parameter should be lower 
than about IQ-^ - 10"^. 

With reactions involving WIMPs in the initial state and standard particles 
as positrons in the final state and supposing these reactions mediated by the U 
boson it would be possible to explain the experimental signals described above. 
Since no astrophysical data involves anomalous production of antiproton, it is 
necessary to require the U boson mass (Mu) to be lower than the mass of two 
protons. 

An interesting consequence of the existence of such a light U boson is that 
it can be directly produced in a controlled environment, such as fixed target 
experiments [15, 16, 18] or high-luminosity e+ e~ colliders at the GeV scale 
(flavor factories) [19-21]. 

At flavor factories, e.g. at DA<I>NE, BESIII and present and future B- 
factories, a particularly clean and simple channel, insensitive to the details of 
how the U boson takes a mass, is surely represented by the resonant radiative 
production of a boson, followed by its decay into a lepton - antilepton pair 
(see Fig. 1). 

A distinctive feature of the expected signal is the appearance of a Breit- 
Wigner peak in the shape of the invariant mass distribution of the lepton pairs 
induced by the mechanism of photon radiative return and corresponding to U 
boson resonant production. The drawback of this channel is the fairly small 
value, over a wide range of parameters, of the signal cross section in comparison 
with the large rate of the backgrounds given by the QED radiative processes 
e+e~ / = e,/i. The signal over background ratio can be enhanced by 
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cutting on the invariant mass of the lepton pair. 

The conclusion presented so far in the hterature is that the production 
process ahows to identify a U boson at present flavor factories if its mass is 
in the range 0.1 — 1 GeV and the kinetic mixing parameter e is greater than 



These analyses are generally based on the evaluation of the number of sig- 
nal events through the calculation of the differential cross section of the 2^2 
process e+e~ Uj^ including the decay of the on-shell U boson into lepton 
pairs by means of branching ratios [19, 22], and/or an approximate estimate 
of the backgrounds. More importantly, all the studies so far performed ne- 
glect the contribution of higher-order initial and final state QED corrections, 
which are known to be a phenomenologically relevant effect at GeV-scale e+e~ 
colliders [23]. 

To investigate whether and how the above approximations affect the physical 
observables of experimental interest, a Monte Carlo event generator, which is 
still missing for such studies, has been made available^ by improving the exist ng 
BabaYaga event generator for QED processes at flavor factories. 

The generator can be useful for data analysis at flavor factories to study 
every physical observable, i.e. invariant mass distributions but also angular 
distributions or correlations and so on. 

It includes: 

• an exact tree-level calculation of the signal and background processes con- 
tributing to the signatures e+e~ ^ 7, /7 ^ /i+/i~7, e+e~7; 

• the effects of the most important higher-order corrections induced by mul- 
tiple photon radiation and vacuum polarization. 

We used our calculation in order to assess its impact on the experimental sen- 
sitivity as evaluated in the literature, and show how this can be enhanced by 
means of event selections not considered so far. 

The calculation has been done by means of the ALPHA [24] algorithm, a 
tool to compute tree- level matrix elements. We implemented the calculation 
in the BabaYaga@NLO event generator [23, 25, 26], to simulate distributions of 
experimental interest and account for realistic event selection criteria. In the U 
boson propagator we include the total dark photon width using the formulae of 
Ref. [20] for the partial widths into leptons and hadrons, i.e. 



where /3f = y 1 — ij/^ R = 1 for leptonic decays of the U boson and 



R = cTe+e- ^hadrons /<^e+e-^/2+/i- fo^" hadrouic decays[27]. It is worth to note 
that the U boson is a dramatically tiny resonance, being is width well under 
any experimental sensitivity, varying from 10~^ to 10~^ for reasonable values 
of Mu and e. 

It is well known that at flavor factories multiple soft and collinear radiation 
emitted by the colliding beams may have a strong impact on the measured 

■"^you can find the code on our website (http://www.pv.infn.it/~hepcomplex). 
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cross section and on the shape of the distributions. The effect of higher-order 
corrections in our Monte Carlo generator is taken into account using the popular 
QED structure function approach [ ]. Initial state radiation modifies the tree- 
level cross section as follows 



dcr(s, t) 



dxi dx2dyidy2daoD'' {xi , s)D'' {x2 , s)D^ {yi , s)D^ (^2 , s) (3) 



where s, i are the Mandelstam invariants after photon radiation and D{x^s) 
is the electron structure function, which describes the probability to find an 
electron with a momentum fraction x inside an electron of a given momentum 
s. If we take into account only the most important terms, due to resummed 
multiple soft and hard photon emission in the collinear approximation 
the structure function has the form: 



D{x,s) 
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il-x)' 



(1 + x) 



(1 + x)(-41n(l - x) + 31nx) -4 



\nx 



(4) 



where P = 2a/7r(ln(5/mj) — 1), F is the Euler gamma function and 7^; is 
the Euler - Mascheroni constant. 

In Eq. 3 D^ixi^ s) refers to the structure functions of the initial state electron 
and D\yi^s) refers to final state fermions. 

Being an effect of size comparable to that of photon radiation, we also con- 
sider in our calculation the running of the electromagnetic coupling constant 
according to (see Ref. [ ] for a recent review) 



Aa{q^) = Aai{q^) ^ Aah{q^) (5) 

where Aai is the contribution due to leptons, which is analitically known, 
and Aah is the non-perturbative hadronic contribution, included according to 
Ref. [2: ] or [-;:]. 

Fig. 2 illustrates a typical signal, the invariant mass distribution of the 
muon pairs for three values of s and Mu at DA<I>NE energies {s = 1.02 GeV), 
as obtained with our calculation in the lowest order approximation. It is possible 
to see a peak in the distribution due to the U boson effects. 



2 Experimental sensitivity 

If we saw a peak in a distribution, for example the one related to the invariant 
mass of final state fermions, we would be able to describe the U boson properties 
like its coupling and its mass. Instead if we couldn't see nothing but the SM 
background, we would exclude a certain region in the Mu - e parameter space. 
So we need to calculate which are the U boson properties that made possible 
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Figure 2: Invariant mass distribution of the muon pairs for different values of 
e and Mu with s = 1.02 GeV 



to confuse the U boson's peak with a statistical fluctuation of the background 
for a certain experiment with its experimental resolution and its luminosity 
{experimental sensitivity). 

On the grounds of the calculation described in Section 1, we revisited the 
experimental sensitivity to a dark force signal evaluated in the literature accord- 
ing to the approximations previously discussed. For concreteness, we consider 
the case of the KLOE/KLOE-2 experiment at the upgraded DA^NE [30] and of 
future experiments as SuperB factories [31, 32]. We also made a comparison with 
BaBar's analysis made by [ ] and the results totally agree, which means that 
other calculations already present in the literature are a really good estimate of 
the exact calculation and that the effects due to the finite width of the U boson 
are negligible for the invariant mass distribution. Neverthless the Monte Carlo 
event generator can be useful for the experimental analysis and for studies of 
backward-forward asymmetry (which is a really smaller signal with respect to 
the presence of the peak in the invariant mass distribution). 

We define the statistical significance {S) as 

5 = ^ = vTi^^-^ (6) 



where (Jf is the full cross section including the exchange of a virtual photon and 
U boson, (Jb {Nb) is the SM cross section (number of events) background, Ns 
is the expected number of events due to the presence of U boson and L is the 
integrated luminosity. In the above equation we require S to be greater than 5 
to claim a discovery. 

We assume L = 5 fb-^ for KLOE/KLOE-2 and L = 100 ab-^ for a B 
SuperB factory, respectively. To simulate detector acceptances, we also impose 
the following energy and angular cuts: 

KLOE 35° < Oi±^^ < 145°, E]f^^ = 10 MeV 
SuperB 30° < 6>^±,^ < 150°, £^[±7^) = 30, (20) MeV 

using as cm. energies ^/s = 1.02, 10.56 GeV, respectively. To improve signal 
events can be detected as peaks in the lepton pair invariant mass close to the 
value M[/, in a window Mu ± Sm- We want Sm to be as small as possible. 
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optimally coinciding with the detector resolution, a crucial parameter for these 
studies. 

As detector resolutions we use Sm = ±1 MeV for KLOE/KLOE-2 and the 
values obtained according to the empirical relations of Ref. [22] for B/SuperB 
factories, giving ^ ±[1 — 10] MeV for a mass Mu in the range 0.1-5 GeV. 

In Fig. 3 we show the reach potential of KLOE/KLOE-2 experiment, where 
the sensitivity to the kinetic mixing parameter e is shown as a function of Mu. 
Actually, we observed that corrections affect with about the same amount both 
the signal and the background cross sections; hence the systematics induced by 
photon radiation largely cancel in the experimental sensitivity and the conclu- 
sions in the presence of QED radiative corrections confirm the ones obtained in 
the lowest order approximation. 




Figure 3: Experimental sensivity for Figure 4: Experimental sensivity for 
KLOE/KLOE-2 SuperB 

Exclusion limits imposed by the precision measurements of the anomalous 
magnetic moment of the leptons and of aQEol^i^^] are indicated as gray areas. 
The most striking features which could be seen in Fig. 3 are the following: 

• the muon channel has a better reach than the e+e~ channel, obviously 
because of the smaller background. The e+e~ channel could instead be 
interesting for low values of Mu, but an experimental difficulty, due to 
the fact that it is really hard to distinguish a low energy electron with 
a photon, exists and limits this kind of analysis. If we want to study a 
U boson lighter than the threshold we have to use some different 
techniques; 

• the sensitivity is significantly degraded if Mu is around the g resonance 
because the branching fraction U is suppressed by the dominant 
decay mode U tt+tt", which would be a more convenient channel in 
this region. 

We have noticed that the maximum sensitivity achievable at the upgraded 
DA^NE with a luminosity L 5 fb"^ , i.e. s - 0.001 - 0.002, is equivalent 
to the one of the present B-factory experiments BaBar/Belle with L 500 
fb~^[-z]. This fact could be easily understood since the reach on e follows the 
rule cx E/^/L. 

In our analysis we realized that radiative corrections do not significatively 
alter the experimental sensitivity because both the background and the signal 
are modified by the same amount. 
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Figure 5: Experimental sensivity for 
KLOE/KLOE-2 - Small Angle selec- 
tion - /i channel 



Figure 6: Experimental sensivity for 
KLOE/KLOE-2 - Small Angle selec- 
tion - e channel 



Fig. 4 presents the experimental sensitivity at a SuperB factory, where it 
could be possible to exclude an epsilon greater than few per 10~^ for all mass 
values up to 2 proton mass. Hence the really high statistics of a SuperB collider 
will allow to probe values of e about an order of magnitude smaller than those 
reachable by present flavor factories. 

Finally, Figs. 5 and 6 show the sensitivity for the small angle selection, 
defined in Eq. 7, for the KL0E/KL0E2 experiment. 

35° < ei± < 145°, EJi''' = lOMeV 

|cos((9^)| > cos(15) (7) 

The small angle selection is also compared with the large angle selection pre- 
viously discussed. It is possible to see that there is a relevant gain in sensitivity 
for Mu greater than ~ 0.5 GeV. 



3 Conclusion 

With existing experiments and analysis of existing data we can exclude a U bo- 
son with an e coupling of about 10~^ for a wide range of mass (flavor factories [ ]) 
and some particular area in the e - Mu parameter space for a really light U boson 
(beam dump experiment [ ]). With a SuperB factory we would be able to see a 
U boson if its kinetic mixing is above few 10~^. Anyway a really big area in the 
e — Mu space would be allowed also after the analysis of a SuperB experiment. 
A U boson with such properties would be able to explain the astrophysical data. 
To cover this area it would be useful a new fixed target experiment [jl i ] in which 
a really high luminosity could be reached. 
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